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Abstract - Anomalous Hall effect arises in systems with both spin-orbit coupling and magne¬ 
tization. Generally, there are three mechanisms contributing to anomalous Hall conductivity: 
intrinsic, side jump, and skew scattering. The standard diagrammatic approach to the anomalous 
Hall effect is limited to computation of ladder diagrams. We demonstrate that this approach 
is insufficient. An important additional contribution comes from diagrams with a single pair of 
intersecting disorder lines. This contribution constitutes an inherent part of skew scattering on 
pairs of closely located defects and essentially modifies previously obtained results for anomalous 
Hall conductivity. We argue that this statement is general and applies to all models of anomalous 
Hall effect. We illustrate it by an explicit calculation for two-dimensional massive Dirac fermions 
with weak disorder. In this case, inclusion of the diagrams with crossed impurity lines reverses 
the sign of the skew scattering term and strongly suppresses the total Hall conductivity at high 
electron concentrations. 


Many ferromagnetic materials exhibit a finite Hall ef¬ 
fect, i.e. transverse voltage in response to a current, with¬ 
out applying external magnetic field. This phenomenon is 
lommonly referred to as the anomalous Hall effect (AHE) 
[l]. Two important ingredients of AHE are magnetiza¬ 
tion and spin-orbit interaction. The former breaks time- 
reversal symmetry and exerts a force acting on electron 
spins while the latter couples the spins to orbital degrees 
rf freedom thus giving rise to the transport effect. 

AHE can also occur as a result of valley or isospin po- 
arization rather than ordinary ferromagnetism [2]. The 
spin-orbit coupling can also be of a more general form as it 
s, e.g. in graphene j[3]|4j where the role of spin is played by 
she sublattice index. An important part of the anomalous 
Hall signal originates in the Berry curvature thus being of 
a topological origin 5|. It is, therefore, natural that the 
discovery of materials like graphene and topological insu¬ 
lators 6,7, which are characterized by non-trivial Berry 
phase of quasiparticles, has considerably broadened the 
interest to AHE from both theory and experiment [§ ■ 14 


Despite the long history 15 17 and high experimental 


relevance of AHE, its theoretical description is a challeng¬ 
ing task often leading to confusions. In modern literature, 


two common approaches based on the Boltzmann kinetic 
equation and Kubo-Streda diagrammatic formalism are 
discussed. Boltzmann equation provides an intuitive qua- 
siclassical approach to the effect 1, 18 j but requires an 
accurate account of several mechanisms of Hall conduc¬ 
tivity: intrinsic, side-jump, and skew-scattering. Intrinsic 
AHE is attributed to topological properties of the band 


19 and is thus independent of disorder. Skew scattering 


is due to the asymmetry in the impurity scattering cross- 
section and side jump refers to the transverse displacement 
of an electron being scattered. An alternative microscopic 
Kubo-Streda formalism is more systematic but less intu¬ 
itive |20[|21|. In this Letter we will employ the latter ap¬ 
proach but also discuss the results in terms of the three 
quasiclassical mechanisms. 


The conductivity of a good metal at low temperature, 
such that direct contributions of inelastic scattering pro¬ 
cesses can be disregarded, can be expanded in powers of 
the large metallic parameter pol (here po is the Fermi mo¬ 
mentum and l is the electron transport mean free path). 
Within the standard model of weak Gaussian disorder, 
the dissipative component of the conductivity tensor be¬ 
haves as oxx oc pol in agreement with the Drude theory, 
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Fig. 1: Diagrams for the anomalous Hall conductivity. Stan¬ 
dard non-crossing approximation (a) yields Eq. CD- The X 
diagram (b) is of the same (zeroth) order in disorder strength 
as all the non-crossing diagrams and corrects the result lead¬ 
ing to Eq. 0- The 'h diagram (c-d) has potentially the same 
order of magnitude but cancels for the massive Dirac model. 
Vertex correction involves the sum of ladder diagrams (e), see 
Eq. {D}. 


while the anomalous Hall conductivity, a xy oc (Po0°> re ~ 
mains constant in the clean limit pol —> oo. Beyond the 
Gaussian approximation, an additional scattering param- 
eter l* < l (related to the third moment of random poten¬ 
tial distribution) arises leading to an enhanced skew scat¬ 
tering, a xy oc PqI* 1.18 22 . In the limit of rare strong 
impurities, this contribution can develop into the resonant 


skew scattering 23 - 26 


In this Letter we consider AHE in the model of weak 
Gaussian disorder using the diagrammatic Kubo-Streda 
formalism 27 . It is well known that in a good metal, 


Pol 3 > h , diagrams with intersecting impurity lines gener- 
ically yield an extra smallness of the order ( pol ) _1 . That 
is why earlier calculations of anomalous Hall conductiv¬ 
ity 20 21 28 employed the non-crossing approximation 


including only the ladder diagram Fig. UK However, we 
show that, contrary to the common knowledge, the “X and 
'L diagrams” with two intersecting impurities, Fig. 0b-d), 
are also crucially important for AHE. Despite impurity- 
line crossing, they contribute to the result in the leading 
order (po0°- 

We argue that importance of these diagrams is a gen¬ 
eral feature of AHE and related phenomena. Indeed, the 
Hall conductivity in the non-crossing approximation con¬ 


tains, in particular, the skew scattering contribution 20 


Since any individual weak impurity has a symmetric scat¬ 
tering profile in the Born approximation, the skew scat¬ 
tering occurs only due to pairs of impurities located at 
a distance less or of the order of the Fermi wave length, 
Po| r — r'| < h. Such closely positioned scatterers cannot 
be treated quasiclassically as isolated objects. Therefore 
a correct treatment of scattering on the double defects 


should necessarily include the X and ’L diagrams. In the 
framework of the Boltzmann kinetic equation, inclusion of 
these crossed diagrams corresponds to a proper calcula¬ 
tion of the full two-impurity scattering amplitude in the 
collision integral. 

To illustrate the general statement formulated above, we 
resort to the minimal model for AHE: 2D massive Dirac 
fermions in the presence of scalar Gaussian disorder, 


H = v erp + ma z + K(r), (la) 

(K(r)K(r')) = 2ira(hv) 2 5(r — r'), (V) = 0. (lb) 

Here er = ( a x , a y ) is the vector of Pauli matrices, p is the 
momentum operator, v is the characteristic velocity, a is a 
dimensionless parameter characterizing disorder strength, 
and the brackets (...) stand for the averaging over disor¬ 
der realizations. 


Aside from random potential, the Hamiltonian (la) con¬ 


tains exactly the two required terms: spin-orbit coupling 
and magnetization. It is thus a generic model of AHE 


28 . Apart from that, such Hamiltonian occurs on the 
surface of a 3D topological insulator subject to a constant 
magnetization m (e.g., due to proximity to a ferromag- 
net). Another possible realization is graphene aligned on 
the surface of BN (in this case to may have large-scale 
spatial variations). The same Hamiltonian describes elec¬ 
tron properties of a narrow HgTe quantum well (linearized 
Bernevig-Hugli-Zhang model 29 ). In the last two cases, 


two mutually time-reversed valleys described by Eq. (la) 
arise, hence cr xy has to be treated as the valley-Hall re¬ 
sponse. For HgTe, it also corresponds to spin-Hall since 
the valleys have opposite magnetic polarization. 

We will evaluate a xy using Kubo-Streda formula 
The Hall conductivity is expressed as a xy = a l x 
with 


xy A ®xy 


<4 = -Tr (a x G R a y G 


A >, 


dn 

= ec 5 B- (2) 


Here G R ' A are the retarded and advanced Green functions 
at the Fermi energy, n is the total electron concentration, 
and B is the external magnetic field. The first term de¬ 
scribes contribution of conduction electrons near the Fermi 
surface while all, accounts for the contribution of the en- 

•by 


tire Fermi sea and has a topological origin 19 . 

When the Fermi energy e is inside the spectral gap, 
|e| < |to|, Hall conductivity is dominated by the second 
term all,, = — (e 2 /2h) signm 5,) 


30 


31 


As the energy is 
increased above to, the conductivity a xy decays from ±1/2 
to zero. We will analyze this dependence as a function of 
the energy and mass (assuming from now on e > m > 0) 
in the leading (zeroth) order in the disorder strength a. 
The topological contribution a l xy is negligible outside the 
gap, aH y c* hence our task is reduced to the analysis of 
the contribution al,,. 

•by 

The standard approach to computing the leading con¬ 
tribution at the Fermi surface is the non-crossing approx¬ 
imation. It amounts to summing up the ladder diagrams 
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Fig. 2: Full anomalous Hall conductivity a xy , Eq. Q, including 
the contribution of the X and 'I' diagrams (solid line). Anoma¬ 
lous Hall conductivity a xy , Eq. (|3|, calculated within the stan¬ 
dard non-crossing approximation (dashed line). 


depicted in Fig. [l](a,e). Solid lines in the diagrams cor¬ 
respond to the disorder-averaged Green function and the 
dashed lines are disorder correlators ©• The logic be¬ 
hind this formalism is based on the fact that the average 
Green function has a sharp maximum at the Fermi mo¬ 
mentum pq = v~ l \J e 2 — to 2 . The contribution of a dia¬ 
gram is maximized when all Green functions are taken at 
p = pq. Each crossing of impurity lines reduces the phase 
space available for the corresponding diagram and thus 
makes it smaller by the factor ( p$l oc a. 

Calculation of a xy for the model (JlJ) in the non-crossing 
approximation was performed in Refs, 
following result: 


20 28 with the 


a xy 


4e 2 em(e 2 


+ m 2 ) 


h (e 2 + 3 to 2 ) 2 


e > m. 


( 3 ) 


This result is plotted in Fig. [2] with the dashed line. In the 
kinetic equation formalism, this Hall conductivity is sepa¬ 
rated into intrinsic, side-jump, and skew-scattering parts 
[see Eqs. (12) and Fig. [4] below]. In the limit e to the 
value of a xy decays as m/e. 

However, as we demonstrate below, the non-crossing ap¬ 
proximation is incomplete. The contribution of X and 'F 
diagrams, Fig.Qb-d), is of the same order a 0 and modifies 
the result leading to 


8e 2 em 3 

"’ = 'T(?TW £>m ’ () 

which is manifestly different from Eq. as illustrated 
in Fig. [2| In particular, the Hall conductivity a xy decays 
as (m/e) 3 in the limit e > m in sharp contrast to the 
behavior of er“. The contribution of X and T diagrams 
is important in the whole parameter range e > m. We 
also show that, for a given model, there exist no other 
contributions to <r xy in the leading order in a. 

Let us now turn to the details of the calculation. We 
start with the computation of the disorder-averaged Green 
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Fig. 3: Electron self energy (a) provides logarithmic renormal¬ 
ization of e and m, while the disorder strength a is renormalized 
due to vertex diagrams (b), see Eq. ©■ Imaginary part of the 
self energy (a) is finite and gives rise to the electron scattering 
rate in the Born approximation. It is explicitly included in the 
average Green function |8]). 


function. For weak disorder, a«l, the lowest order self 
energy (Born approximation) is given by the diagram de¬ 
picted in Fig. UK The real part of this diagram diverges 
logarithmically with the band width Ao- This divergence 
indicates that the Born approximation may be insufficient 
for the massive Dirac model ([l]). Instead, we perform an 
accurate summation of leading logarithmic contributions 
with the help of the renormalization group (RG) proce¬ 
dure [30 32-34 . Apart from the self energy diagram Fig. 


UK RG also includes logarithmically divergent vertex di¬ 
agrams Fig. UK As a result of renormalization, the pa¬ 
rameters of the model depend on the value of the running 
ultraviolet cutoff A, 


da „ , dine din to. 

= 2 a, —— = at, — r = -a. (5) 


din A 


din A 


din A 


Starting from the band width Ao and bare parameters of 
the model, which are temporarily denoted as eo, Too, and 
ao, we solve the one-loop RG equations ([5]) and obtain the 
observable (renormalized) values 

a = Za o, e = Z l ^ 2 eo, to. = Z~ 1 ^ 2 mo, (6) 

Z(A) = [l-2a 0 ln(A 0 /A)]- 1 . (7) 

The RG process has to stop when the running scale reaches 
A = V e 2 — to 2 . At this point all logarithmic contribu¬ 
tions are taken into account. We assume that the disorder 
strength a is still small after the renormalization. 

From now on we use only renormalized parameters of 
the model and hence disregard the real part of the Born 
self energy Fig. UK The imaginary part of the self en¬ 
ergy results in small imaginary contributions to energy 
and mass yielding the retarded Green function 

R _ e(l + iira/2) + ma z ( 1 - iira/2) + erp 
Gp ~~ e 2 — to 2 — p 2 + ina(e 2 + to 2 ) ' W 

Advanced Green function is the hermitian conjugate of 
the above expression. Here and below we let H = v = 1 to 
simplify notations. 

We use the averaged Green function ([8]) to calculate a 
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single-impurity vertex correction, 
72 „ 


2na j Jfrr^ G v axG v = A(Tx + BrJy ’ 


A = 


2(e 2 


B = - 




m 


2 ’ 


(9) 

( 10 ) 


Note that the diagonal term A is of the zeroth order in the 
parameter a while the factor B oc a is small. Summing up 
the vertex ladder, see Fig. [ij, we include the B term only 
once. This yields the following expression for the dressed 
current operator: 


a x = Fa x + F 2 B<j v , F = 


1 


1 - A 


2(e 2 + to 2 ) 
e 2 + 3 to 2 


( 11 ) 


The non-crossing approximation to a xy is the result of 
adding one last ladder rung § to the already dressed 
current vertex d x and taking the trace with the bare cur¬ 
rent operator a y . This results in a xy = F 2 B/ira, which 
coincides with Eq. (J3| obtained previously in Ref. 28 


To understand why the result of Eq. <© is incomplete, 
one has to figure out which characteristic momenta pro¬ 
vide the main contribution to it. This has been already 
investigated in detail in Refs. 18 20 , where the same 


problem had been addressed in the eigenbasis of the bare 
Hamiltonian. This basis contains two eigenstates |±) per 
each momentum corresponding to the two bands of the 
Dirac Hamiltonian ± \Jp 2 + m 2 . Velocity operator does 
not commute with the Hamiltonian and hence has finite 
off-diagonal matrix elements in this basis. Setting all the 
Green functions on the mass shell e = \Jp 2 A to 2 (i.e., 
completely disregarding the negative “positron” band) 
correctly reproduces cr xx but leads to vanishing a xy . 

In order to obtain a finite Hall conductivity in the ze¬ 
roth order in a, one has to allow for exactly one internal 
momentum in the diagram laying outside the mass shell. 
In the bare bubble with no impurities, this is achieved 
by taking one Green function, which connects the two 
current vertices, away from the mass shell (i.e., two off- 
diagonal matrix elements of velocity are involved). This 
yields the intrinsic contribution to AHE. Similarly, the 
side-jump term corresponds to diagrams containing only 
one off-diagonal matrix element of velocity. In these dia¬ 
grams, the off-mass-shell Green function connects the off- 
diagonal current vertex to an impurity. Finally, the skew¬ 
scattering diagrams contain two diagonal current vertices 
while the off-mass-shell Green function connects two im¬ 
purities. Within the non-crossing approximation, these 
three components of the anomalous Hall conductivity are 


given by 20 28 


int_ 

a xy ~ 


(T SJ = 
xy 


ss-nc _ 

a xy ~ 


h 2e ’ 

e 2 2to(c 2 — to 2 ) 
h e(e 2 + 3 to 2 ) 
e 2 3 m(e 2 — to 2 ) 2 
h 2e(e 2 + 3to 2 ) 2 


(12a) 

(12b) 

(12c) 



Fig. 4: Intrinsic (12aI, side-jump (12bI, and skew-scattering 
(181 components of the anomalous Hall conductivity (solid 
lines). Dashed line shows the skew scattering part in the non¬ 
crossing approximation (12c I. 


They are plotted in Fig. [4j The Green function taken 
away from the mass shell decays rapidly in real space. 
Therefore the latter skew-scattering contribution is due to 
close pairs of impurities (separated by a distance d < p^ 1 
smaller or of the order of the Fermi wave length). This 
means that crossing the corresponding impurity lines does 
not introduce an extra smallness! That is why the X and 'F 
diagrams depicted in Fig. [ljb-d) are of the same (zeroth) 
order in a as all the non-crossing terms and provide an 
important additional part of skew scattering. 

In order to calculate the contribution of X and 4/ di¬ 
agrams we resort to the real space representation. The 
diagrams of Fig. 0b -d) are given by the following inte¬ 
grals: 

4 = ^ J d 2 rTr( y J*G B _ r jyG^ r ' S ), (13a) 

a* y = ^ J d 2 r T \(j*G B r G B J v _ r A h.c.), (13b) 

where the vector r connects the impurities. We define J r 
as the Fourier transform of the dressed current operator 
together with two adjacent Green functions: 


T =F [*L- C a <tC b F pr = 2V [^-^] 
r F J (2ir) 2 Gp Gp 7ra(e 2 + 3 m 2 ) 


(14) 


The last expression holds to the leading order in a. The 
two elements containing J in each diagram yield aT 2 that 
cancels the contribution of the two crossing impurity lines. 
Hence upon substitution of Eq. (14) into Eqs. (13), all 


Green functions can be considered in the clean limit. This 
also validates the use of the diagonal part F of the dressed 
current operator in Eq. ( fl4| . 

The clean Green function has the following form in real 
space: 


G B ’ A = 


- [e + ma z - icrV] |Yo(P(R) T iJoipor)]. 


(15) 


where Jq and Yo are the Bessel functions of the first and 
second kind, respectively, and po = Ve 2 — m 2 is the Fermi 
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momentum. After taking the trace and integrating over 


directions of r, the integrals in Eqs. (13) reduce to 


that case, X and \F diagrams lead to logarithmic renor¬ 
malization : 30,32-341 of a but are excluded from the weak 


x 47re 2 em(e 2 — m 2 ) 

a xy 


h (e 2 + 3?n 2 ) 2 
27re 2 em(e 2 — m 2 ) 
xy ~ Hr (e 2 + 3m 2 ) 2 


al, = 


JdsJfiJiYo- J 0 Yi), (16a) 

J ds [JUJiYo + J 0 Yi) 

+ J 0 Jiy 0 (Jo- j 2 )], (16b) 


where the argument s = por of the Bessel functions is sup¬ 
pressed. Note that the main contribution to the integrals 
in Eqs. (161 is indeed provided by s ~ 1, i.e., by distances 


localization since they do not contribute to magnetoresis¬ 
tance (do not form a loop in real space). Our result for 
the Hall conductivity Eq. © does not explicitly depend on 
a hence the renormalization of disorder strength is unim¬ 
portant. Computing weak localization correction to a xy , 
however, will require inclusion of additional X- and 'F-type 
diagrams in the Cooperon loop. This will be the subject 
of a separate publication. 

The massive Dirac Hamiltonian (la) possesses non- 


of the order of the Fermi wave length l/po as anticipated. 
Evaluation of Eqs. (16) yields the result: 


4e 2 


i(e 2 — m 2 ) 


'xy 


h (e 2 + 3m 2 ) 


2 ’ 


aZy = 0 . 


(17) 


trivial topological properties that are manifested in the 
finite Berry phase 7r accumulated along a large loop in 
momentum space. This is the origin of many fascinating 
topological effects including odd quantization of the Hall 
conductivity in graphene [3j[4], quantum spin-Hall effect 
in HgTe [8], and topological magnetoelectric response at 


We see that the T diagram actually vanishes while the X 
term gives a positive contribution in the skew-scattering 
channel. Cancellation of the *F diagram is accidental for 
the massive Dirac model and does not hold for a more gen¬ 
eral model of AHE. Comparing to the results of Ref. 28 


the surface of topological insulators |6,7 37. It is also 


we conclude that the whole skew-scattering part changes 
sign upon inclusion of the X diagram (see Fig. El), 


e 2 (e 2 — m 2 )(5e 2 + 3m 2 ) 
h 2e(e 2 + 3m 2 ) 2 


(18) 


The full result for the anomalous Hall conductivity in the 
zeroth order in a is given by Eq. Q. 

Thus we have shown that the X diagram contributes to 
the anomalous Hall effect in the leading order with respect 
to the disorder strength. (This result has been anticipated 
already in Ref. 


intimately related to the fermion doubling. A finite Berry 
phase requires the presence of the second valley in the 
spectrum or a second (opposite) surface in the case of the 
topological insulator. The Hall conductivity, observable in 
a transport experiment, will necessarily include the con¬ 
tribution of both valleys or surfaces 31 36 . Each of these 
contributions will have the form of Eq. (|4j) possibly with 
different values of e and m. 

Other models of AHE include kinetic (quadratic in mo¬ 
mentum) terms in the Hamiltonian. Two most studied 
models of this type are the ferromagnetic Rashba model 
[l] and the Bernevig-Hughes-Zliang model 29 . Bending 


of the spectrum due to quadratic terms removes the non¬ 
trivial total Berry phase but doubles the Fermi surface. 
This significantly complicates the calculation of the Hall 


20 but the contribution was incorrectly conductivity 21 as compared to the massive Dirac model 


regarded as being parametrically small.) While evaluation 
of the skew scattering on double-impurity requires an ac¬ 
curate quantum-mechanical calculation including X and 'F 
diagrams, the resulting cross-section can be directly used 
in the collision integral within the quasiclassical Boltz¬ 
mann kinetic equation. The above result is then repro¬ 
duced along the lines of Ref. 


considered here. Nevertheless, the X and diagrams must 
be also included in a xy within these general models. The 
result of this calculation will be published elsewhere. 

Skew scattering is relevant for a number of related phe¬ 
nomena such as spin-Hall effect [2], polar Kerr effect in 
chiral superconductors 38-40 , spin-orbit torque 41 42 


20 


When three or more intersecting impurity lines are in¬ 
cluded in the conductivity diagram, they produce an ex¬ 
tra smallness leading to the weak localization correction 
Sa xy oc a. The general logic in estimating the order of 
a diagram is the following. Finite Hall conductivity re¬ 
quires at least one Green function in the diagram to be 
taken away from the mass shell p = po- The impurities 
connected by this Green function have to be considered as 
a single scattering complex. Once these two impurities are 
combined, counting remaining crossings determines the or¬ 
der of the diagram in the parameter a. Equivalently, the 
order of the diagram can be estimated from the number of 
loops in real space after the two close impurities are fused 
together. 

Weak localization correction to the longitudinal conduc¬ 
tivity a xx for the model ([l]) was computed in Ref. 35 . In 


etc. A proper quantitative description of such effects also 
requires inclusion of X and ^F diagrams. 

In conclusion, we have demonstrated that an accurate 
evaluation of the anomalous Hall conductivity should in¬ 
clude the diagrams with two crossed impurity lines, Fig. 
[ljb-d), in addition to the standard set of non-crossing di¬ 
agrams, Fig. EH Despite impurity-line crossing, these 
diagrams contribute in the leading order to the skew scat¬ 
tering component of anomalous Hall conductivity. In the 
model of weak disorder, skew scattering originates solely 
from pairs of impurities at a distance smaller or of the or¬ 
der of the electron wavelength. Such double defects require 
full quantum-mechanical treatment beyond the commonly 
employed ladder approximation. We illustrate our findings 
using AHE with two-dimensional massive Dirac fermions 
as a model example. Upon inclusion of X and 4^ diagrams, 
the anomalous Hall conductivity © is considerably sup- 
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pressed as compared to the non-crossing approximation 
result ([3]). We argue that X and T diagrams are indis¬ 
pensable for all models of AHE as well as for a number of 
related phenomena. 
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